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ABSTRACT 

We used Optical, Near Infrared photometry and radial velocity data for a sample of 11 
Cepheids belonging to the young LMC blue populous cluster NGC 1866 to estimate their radii 
and distances on the basis of the CORS Baade-Wesselink method. This technique, based on an 
accurate calibration of the surface brightness as a function of (U-B), (V-K) colors, allows us to 
estimate, simultaneously, the linear radius and the angular diameter of Cepheid variables, and 
consequently to derive their distance. A rigorous error estimate on radius and distances was 
derived by using Monte Carlo simulations. Our analysis gives a distance modulus for NGC 1866 
of 18.51±0.03 mag, which is in agreement with several independent results. 

Subject headings: star clusters: individual(NGC 1866) — variable stars: Cepheids — distance scale 
Introduction 



accuracy of 1-2% (jFreedman et al.ll201ll . and ref- 
One of the main goal of modern cosmology is erences therein). The distance to the LMC is a 

the determination of the Hubble constant to an critical step in the problem of determining the 



1 



scale of the universe. It is, in fact, considered as 
a benchmark in the determination of distances to 
other galaxies, being a useful place to compare 
and test different distance estimators. Therefore, 
any error in its distance contributes a substantial 
fraction of the unce rtainty in the Hubble constant 
(|Mould et al.ll200fl) . 

In the last two decades there have been im- 
provements in the determination of the LMC dis- 
tance both theoretically (e.g. red clump, tip of red 
giant b ranch, model fitting of variable stars light 
curves; Romaniello et al.ll2000 : Walker et al.ll200lt 



Keller fc Woodll2002l:lBono. Castellani fc Marcoi 
20021 : iMarconi fe ClementinJl2005tlKeller fe Woo( 



20061; iGrocholski et all 120071; iKoerwerl 
and observationally (e.g. parallaxes for Galac- 



PoDOwski & Gouldlll998: Whitelock fcFeastl 


2000; 


Benedict et al.l2002t FitzDatrick et alJ 20021 


2003 


DalPOra et al.l 12004 Benedict et al.l 120071 


2011 


Bonanos et al.l 20111). The literature contains a 



huge number of esti mated distanc e values rang- 
i ng from 18.10 mag (IUdalskilll998l) to 18.80 mag 
([Groenewegen fc Oudmaijerl |2000( V a very unsat- 
isfactory situation that indicates the presence of 
significant systematic errors in most of the meth- 
ods. A summary of the results of different methods 
with the correspondin g references can be found in 
Benedict et all ([20021 Tab. 10). 



Since the discovery bv lLeavitt (1908), the most 



used absolute calibration of the extragalactic dis- 
tance scale is based on the Period-Luminosity 
relation of LMC Cepheids. Unfortunately, this 
procedure typically gives an indirect measure of 
the distance, because it calibrates the Period- 
Luminosity relation for LMC on the basis of the 
sam e relation obta i ned fo r Galactic Cepheids (see 



e.g. IPersson et al. 2004 ). A strong drawback 



in this procedure is the role played by mctal- 
licity in the Period-Luminosity relation. From 
both the observational and theoretical perspec- 
tive, some au thors argue that this dependence is 
prese n t (e.g. iFreedman et al.l l200lfc IStorm et al 

' 20081: 



2003; Tammann Sandaec fc Reindl et al 



2003; Gieren et al 



Alibert et al.lll999h 



20051 : IStorm et al.ll2005l [201 lbt 



In the Hubble Space Te lescope H Key Project 
by IFreedman et all ( 2001 ) a LMC distance mod- 
ulus value of 18.50±0.10 mag has been adopted 
to calibrate several secondary distance indicators 
using the Cepheid variables. The adopted value 
of the LMC distance by Freedman and collabora- 
tors can be considered as a converging value for 
the LMC distance and as the watershed for the 
dichotomy between the so called "short distance 
scale" (values lower than 18.50 mag) and the "long 
distance scale" (values larger than 18.50 mag). 

A direct measurement of the distance to the 
LMC can be obtained by using a version of 
the Baade-Wesselink te chnique called the sur- 



face brightness method (jGieren et al.l 119971 and 
references therein). This method utilizes radial 
changes of the stellar surface brightness and the 
pulsational velocity for the determination of the 
linear r adius and angular d iameter of pulsating 



stars. Gieren et al 



(|2005l ) applied the near- 
infrared surface brightness technique to a sam- 
ple of 13 Cepheids in the LMC and obtained a 
distance modulus of 18.56±0.04 ma g. The same 
method has been applied recently by IStorm et al 



(|2011bt ) who derived the distances of 36 Cepheids 
in the LMC and obtained a distance modulus of 
18.45±0.04 mag. 

In the present work we are going to face 
the problem of the LMC distance by means of 
the CORS version of the Baade-Wesselink sur- 
face brightn ess tec hnique 



Ripepi et al.1 



Molinaro et al 



9971 hood. 



(Caccin et al. 
IRuoppo et al.l 



1981 



2011) 



2004; 



using recent photometric 
and spectroscopic data of a sample of Cepheids 
observed in the populous cluster NGC 1866. 
Th is object is one of the few young (~ 10 8 



Brocato et al. 19891 ) LMC clusters that are 



yr 

close enough to allow the det ailed observation 



of individual stars (see e.g. iMucciarelli et al.l 



2011 : Brocato et al. 20031 an d references therein) 



Moreover. lStorm et al.1 (|2005l ) found that the NGC 
1866 Cepheids are close to the average LMC 



Sakai et al.ll2004l:lRomaniello et al.l2005l:lMacri et al.l r^pi^ri distance, so we will assume the distance 



2006 : lFreedmanfc Madord201ll:lMarconi. Musella fc 



20051; iBono et alll2008l; iBono et al.ll2010h , even if 
the size and the sign of the effect is still disputed, 
while others suggest that the effect is small and 
perhaps ill defined (e.g. Groenewegen fc Salarid 



Fior^tiflsls cluster to be the distance of the LMC 
itself. To date it is known that NGC 1866 har- 



bors at least 23 Cep heids (|Welch fc Stetsonlll993 : 



Musella et al.l 120061 ). the largest number among 



all the LMC clusters, and this makes it an excel- 



2 



lent system for distance determination. Indeed it 
provides a sizeable sample of variable stars con- 
tained in a limited volume and then all at the 
same distance. Furthermore, we expect they are 
characterized by the same chemical composition, 
age and reddening, so that we can derive their dis- 
tance with no influence of differences in the above 
quantities among stars. 

The first Baade- Wesselink distanc e to NGC 
1866 was derived by ICote et al. ( 1991 ). obtaining 
the radii of seven variables in the field of NGC 
1866 and using the B-V color to calibrate the sur- 
face brightness. The final result, 18.6±0.3 mag, is 
affected by a large error, probably due to the used 
color index , which is not a good surface brigh tness 
calibrator ( Coulson Caldwell fc Gieren 19861) . 

Gieren. Richtler fc Hilkerl (|l994l ). using the 



sions are contained in Sec. [7] 



2. The data 



surface brightness modification of the classical 
Baade-Wesselink method, obtained the distances 
of four Cepheids in the field of NGC 1866. They 
used the V-R color index and the radial velocities 
obtained from spectroscopic data tak en at the Las 



Campanas 2.5 m du Pont reflector ([Welch et al 
1991). The obtained distance modulus to NGC 



1866 was 18.47±0.20 mag, considering only three 
Cepheids (HV 12198, HV 12199 and HV 12203). 
The fourth, HV 12204, was excluded because the 
authors suspected it was not a member of the 
cluster. 

Finally, using o ptical and near-infrared data, 
Storm et all (|2005l ) derived the distance to NGC 
1866 through the surface brightness method ap- 
plied on five Cepheids and found 18.30±0.05 mag. 
We aim at improving their result by increasing the 
analyzed sample of Cepheids, thanks to new pho- 
tometric and spectroscopic data, and by relying on 
an accurate calibration of the surface brightness 
obtained from grids of theoretical atmospheres. 

Sec. [2] contains a description of the photomet- 
ric and spectroscopic data used in this work. The 
procedures followed to phase the light curve the 
radial velocity curve and to correct for redden- 
ing are described in Sec. El The CORS Baade- 
Wesselink method is introduced in Sec. |4] together 
with the procedure adopted to calibrate the sur- 
face brightness function using atmosphere models. 
Sec. [5] contains the derived Cepheid angular diam- 
eter and the linear radii. The distance to NGC 
1866 and the comparison with other results from 
the literature are discussed in Sec. [6] while conclu- 



N GC 1866 harbors at le ast 23 Cepheids (|Welch k, Stetson 
19931 iMusella et alj|2003) . Among them we have 
selected those stars such that both near infrared 
and radial velocity data were available. The final 
sample consists of 11 Cepheids whose data are pre- 
sented in the next two sections and are resumed 
in Tab. HJ 

2.1. Photometry 

We used photometric data in the ultraviolet U 
band, optical B, V, I bands and near infrared K 
band. 

The U band data were obtained at Cerro Tololo 
Inter-American Observatory (CTIO) with 0.9 m, 
1.5 m and 4.0 m telescopes. 

The optical data were taken at VLT with 
FORS1 imager and consist in 69 images in B, 
90 images in V and 62 images in I. Both ultra- 
violet and optical data were calibrated in the 
Johnson-Cousins photometric syste m. These VLT 
data h ave already been published in lMusella et al.l 
(|2006l) . The complete UBVRI dataset used in this 
work, will be published in a forthcoming paper 
(Musella et al. in preparation), including a de- 
tailed discussion of the adopted procedure to re- 
duce the data and to face the crowding problem. 

As for the near infrared K band, we have used 



the dataset descri bed in IStorm et al] (|2005) and 
Testa et al.l ([20071 ). Furthermore, we used unpub- 
lished data collected by J. Storm for the Cepheids 
V4, V7 and V8, wh i ch we re calibrated using the 
data by Testa et al.l ( 2007 ) as reference. The data 
from iTesta et al.l (|2007t) were calibrated in the 
LCO phot o metri c system using the r elations from 
Carpenter ( 2001 ). K band data from lStorm et all 
(2005) lave been calibrated in the CIT photomet- 
ric system. 

To compare the photometric data with theoret- 
ical models (see Sec. l4.2|) we have transformed the 
infrared data from CIT and LCO to SAAO photo- 
metr ic system through the simple transformation 
from iBessell fc Bret? (jl988l) : 



Ksaao = Kcit + 0.014 = K LCO + 0.014 (1) 

where we have used also the relation Kcit = 
Klco from lStorm et alJ (|2005h . 
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2.2. Radial velocity 



We used the radial veloc ity d ata fro m Storm et al 
(|2005l ). IStorm et all i2004 ) and lWelch et al.l (Il99ll) . 
obtained through the classical cross-correlation 
method. For three Cepheids, namely HV 12197, 
HV 12199 and We2, we used ne w data coming 
from t he FLAMES VLT dataset bv lMucciarelli et al 
(l201lL Appendix©. We note that our radial ve- 
locities are the first ones ever published for We2. 
The radial velocities for the quoted three Cepheids 
are listed in Tab. [TJ 

The new radial velocity curves for HV 12197, 
HV 12199 and We2 are shown in Fig. [1] together 
with those from the other authors. The agreement 
between our data and the other samples is evident 
for HV 12197 and HV 12199, while a well-covered 
radial velocity curve is shown for We2. As for We2, 
we have also estimated the systematic velocity V 1 
obtaining 301.4 km/s with an accuracy ~1 km/s, 
confirming its cluster membership. 

Finally we recall that from these spectroscopic 
data the iron content was also estimated following 



the approach described in lMucciarelli et al.1 (|2011l ) 
(see also Appendix |A"| . The [Fe/H] values for the 
three Cepheids are -0.39±0.05 dex for HV 12197, 
-0.38±0.06 dex for HV 12199 and -0.43±0.05 dex 
for We2. These values are fully consistent with 
the average iron content of the cluster [Fe/H]=- 
0.43±0.01 dex obtained from the analysis of a sam- 
ple o f static stars in NGC 1866 ( Mucciarelli et al.l 
201 lh . 



H — ' — ' — 1 — 1 — h 

HV 12199 



H — 1 — 1 — 1 — 1 — h 



Fig. 1. — The new radial velocity data for HV 
12197, HV 12199 and W e2 (filled circles) are shown 
together with those from Storm et al.l (|2005l ) (open 
triangles) and Welch et al.l ( 199lh (open squares). 
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Table 1: The new radial velocities in Km/s for HV 12197, HV 12199 and We2 are listed in second, third and 
fouth columns respectively. The epochs of observations arc in the first column. 



Epoch (JD) HV 12197 HV 12199 We2 



53280 


.77407 


312. 


.8±0. 


.8 


320 


,3±1 


.0 


276 


1± 


1 


.2 


53335 


.71398 


287. 


.8±0. 


.7 


314. 


,8±1. 


.1 


285 


,6± 


1 


.0 


53335 


.75731 


289. 


.8±1. 


.1 


315. 


,0±0. 


.8 


285 


,0± 


1. 


.0 


53335 


.80046 


289. 


.8±0. 


.6 


316 


,8±0. 


.9 


282 


,6± 


0. 


.9 


53338 


.60509 


283. 


0±0. 


.7 


319. 


,8±0. 


.9 


302 


.2± 


1. 


.0 


53351 


.59229 


293 


.5±1. 


.0 


316 


,9±1. 


.0 


286 


,4± 


1 


.1 


53351 


.64296 


293. 


.4±1. 


.1 


318. 


,1±1. 


.2 


286. 


,3± 


0. 


.8 


53351 


.69311 


294. 


1±0. 


.9 


319. 


,4±1. 


,0 


288 


1± 


0. 


.9 


53376 


.61793 


290. 


.8±0. 


.8 


282. 


,1±0 


.8 


301 


.4± 


1. 





53376 


.72984 


292. 


.7±1. 


.0 


284. 


,4±0. 


.9 


302 


,9± 


0. 


.8 


53377 


.73544 


312. 


.7±1. 


.1 


311. 


,7±1. 


.0 


321 


,4± 


0. 


9 


53378 


.59190 


298. 


1±0. 


.8 


312 


,3±1 


,0 


282 


,1± 


0. 


.8 
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Table 2: The Cepheids considered in this work are listed in the first column, their period is in the second column and in the third column 
we list the source of the radial velocity data with the shift (Km/s), if any, in parentheses. Finally the epochs used to phase light curves 
and radial velocity curves are listed in the last column. 



Star 


Period (days 


) Radial velocity references 


Epochs (HJD) 


HV 12197 


3.143742 


Storm et al. (2005). this workf-1.251. Welch et al. (1991) 


2452998.5638 


HV 12198 


3.522805 


Storm et al. (2004). Storm et al. (2005) (+2.8). Welch et al. (1991)(+1.1) 


2452947.7500 


HV 12199 


2.639181 


Storm et al. (2005). Welch et al. f 1991V this work 


2452993.6100 


HV 12202 


3.101207 


Storm et al. (2005). Welch et al. (1991) a 


2452930.8922 


HV 12203 


2.954104 


Storm et al. (2005). Welch et al. (1991) 


2452972.6731 


V4 


3.31808 


Storm et al. (2005). Welch et al. (1991) 


2452998.5651 


V6 


1.944252 


Storm et al. (2005). Welch et al. (1991) 


2452976.7400 


V7 


3.452075 


Storm et al. (2005)(-3.4). Welch et al. (1991) 


2452915.9000 


V8 


2.007157 


Storm et al. (2005). Welch et al. (1991) 


2452915.7700 


We2 


3.054847 


this work 


2452990.6712 


We8 


3.039855 


Welch et al. (1991) 


2452995.5394 


"The shifts arc 


the same as those described inlStorm et al.l l|2005f). 





3. Data analysis 

In this section we describe the steps performed 
to phase the light and the radial velocity curves, as 
well as to estimate the reddening in the direction 
of NGC 1866. 

3.1. Construction of light curves and ra- 
dial velocity curves 

Cepheid photometric data were phased as 
usual, i.e. requiring that the maximum of light 
in B band occurs at phase zero, while the max- 
imum in t he other bands are shifted in phase as 
expected (lLabhardt. Sandage fc Tammaim 1997 ; 
iFreedmanl Il988l see e.g.). The adopted epochs 
and the periods are listed in Tab. [2j As the data 
were obtained at different times, to estimate the 
(V-K) and (U-B) colors, it was necessary to first 
interpolate each light curve at the same phase 
points. This was achieved by means of a C code, 
written by one of the authors, performing smooth- 
ing spline interpolation. A different procedure has 
been performed to interpolate the U band of HV 
12197 and K band of We2 and We8. HV 12197 
U band light curve, in fact, was poorly sampled, 
while the K band measurements for We2 and We8 
were of lower quality with respect to the other 
variables. In these three cases we have performed 
the interpolation by using a template light curve. 
To construct the template we have calculated the 
mean ratios A(K)/A(V) and A(U)/A(V) of the 
light curve amplitudes in the K, U and V bands for 
all the stars of our sample showing well-covered 
light curves. Considering the V band light curve 
of the star to be interpolated, we have then cal- 
culated a "normalized" template light curve char- 
acterized by zero mean magnitude and amplitude 
equal to one. Finally, the interpolation has been 
achieved by rescaling the template light curve, ac- 
cording to the estimated ratios, and shifting it to 
the mean magnitude of the investigated Cepheid. 

To estimate the uncertainties of the interpo- 
lated magnitudes for all the Cepheids we have con- 
sidered the rms of the residuals around the inter- 
polated curve. 

As for radial velocity curves, to obtain an accu- 
rate phasing with light curves, we have used the 
same epoch and period as for the photometry. 

As cited in Sec. 12.21 for all the Cepheids, ex- 
cept We8, we have measurements from different 



data sets. To combine them, we corrected for pos- 
sible shifts in radial velocity (see Tab. [2]). They 
were determined by interpolating the most accu- 
rate sample of radial velocities and then by es- 
timating the median distance along the velocity 
axis between the fitted curve and the other radial 
velocity samples. 

3.2. Reddening determination 

To correct the colors for extinction , we used the 



grids of models bv lBessell et al.l ([1998') (http://wwwuser.oat.ts.astro.it 



They provided models with different metallicity 
values and we have linearly interpolated between 
solar ([Fe/H]=0.0 dex) and subsolar ([Fe/H]=- 
0.5 dex) metallicity grids to obtain models at the 
metallicity estimated by iMucciarelli et al.l (|201ll) 
and equal to [Fc/H]=-0.43 dex. 

In this step we fol l owed the technique intro- 



duced by iDean et al 



(|l978l) to derive the red- 
dening of Cepheids. According to this method, 
the reddening is derived by shifting the observed 
points in a color-color diagram onto their intrin- 
sic unreddened position, which can be obtained 
either using stars of known reddening or by using 
models. In particular, according to the models, 
the Cepheids occupy a narrow zone in the plane 
(V-I)-(B-V) (Fig. [2J), which represents the unred- 
dened locus (5000< T e <6250, 1.0< \ogg <4.0). 
Using this color-color plane, we plotted the point, 
whose coordinate are the mean observed colors of 
each Cepheid, and shifted it along the reddening 
vector, trying different values of E(B-V) and se- 
lecting all the values which make the star to lie on 
the grids. The mean of these values is an estimate 
of the E(B-V), and the width of the range of se- 
lected values can be assumed as the uncertainty on 
the mean value. The procedure is shown in Fig. [2] 
for the Cepheid V7 and the values of E(B-V) for 
all the other stars are listed in Tab. [3] We re- 
peated the procedure using t he unreddened locus 
defined by IDean et al. ( 1978f) finding an excellent 
agreement. 

The procedure does not work for HV 12197 be- 
cause its representative point is bluer than the 
grids along the V-I color. A possible explanation 
of this behavior can be the presence of a compan- 
ion which can affect the color of the Cep heid expe- 
cially in optical bands ( Szabadosl 2003 ). In order 
to reduce the effect of outliers, we used the median 
of the derived reddening values obtaining E(B- 
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Table 3: The values of reddening E(B-V) (second column) and their uncertainties (third column) for the 
selected Cepheids of NGC 1866 (first column). 



Star E(B-V) <5E(B-V) 



HV 12197 






HV 12198 


0.102 


0.018 


HV 12199 


0.082 


0.018 


HV 12202 


0.03 


0.02 


HV 12203 


0.02 


0.02 


V4 


0.07 


0.02 


V6 


0.05 


0.02 


V7 


0.06 


0.02 


V8 


0.091 


0.009 


We2 


0.10 


0.03 


We8 


0.06 


0.02 




Fig. 2. — A schematic plot of the procedure fol- 
lowed to estimate the reddening E(B-V) from grids 
of theoretical models in the (B-V)-(V-I) plane. 
The locus expected to be occupied by Cepheids 
is enclosed between constant gravity and constant 
temperature curves (black and red lines respec- 
tively). As example, the representative point of 
the Cepheid V7 is out of the grids when it is not 
corrected for reddening. The values of E(B-V) 
which make the point to move along the reddening 
vector on the grids are selected (dashed line). 



V)=0.06 mag and a rms of 0.02 mag, consistent 
wit h the value typica lly adopted for NGC 1866 
(see IStorm et al. 20051 and references therein). 

To calculate the extinction in each observa- 
tional band, we used E (B-V)=0.06 m ag, Ry = 3.3 
Feast fc Walkerll987t ) and the law bv lCardelli et al 
1989J) obtaining A v = 0.303 mag, A B = 0.261 
mag, Ay = 0.198 mag, Aj = 0.120 mag and 
A K = 0.02 mag. 

4. The CORS Baade Wesselink method 

In the following sections, we give a brief de- 
scription of the complete CORS Baade-Wesselink 
method used to derive the linear radius of Cepheid 
stars and of the procedure followed to calibrate the 
surface brightness function, which is at the base 
of the complete CORS technique. Moreover, the 
surface brightness allows us to estimate the star 
angular diameter and consequently the distance d 
which is the factor linking the linear radius R and 
the angular diameter 9 according to R = \dd. 

4.1. Theoretical background 



Xhe original CORS method (jCaccin et al 
is a realization of the classical Baade- 



1981) 



Wesselink method ( Wesselink 19461) useful to de- 
rive the radius of pulsating stars. 

Starting from the surface brightness function: 

S v = my + 5\og8 (2) 

where 9 is the angular diameter (in mas) of the 
star and my is the apparent visual magnitude, it 
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is straightforward to obtain the basic CORS equa- 
tion by differentiating with respect to the phase 
(</>) , multiplying the result by a generic color index 
(Cij) and integrating along the pulsational cycle: 



-B + AB = 
( 3 ) 

where q — jj^q, P is the period, v is the radial ve- 
locity and p is the radial velocity projection factor, 
which correlates radial and pulsational velocities 
according to R'((j>) — —p-P- v((f>). The last two 
terms, B and AB, are: 



q [ ]n\Ro-pP [ v(<f>')d<f>'}c' i:j d<f> 

JO Jcf>o ' 



B 



Cij{(f))m' v ((f))d(f) 



AB = / Ci^cfyS'yi^dcj) 



(4) 
(5) 



Equation ([3]) is an implicit equation in the un- 
known radius Rq at an arbitrary phase </> . The 
radius at any phase <f> can be obtained by inte- 
grating the radial velocity curve. The main char- 
acteristic in the Eq. © is the estimate of the AB 
term as it contains the surface brightness func- 
tion. Typically the AB term has a small value 

(io- 3 



10 4 , lOnnembo et al.l 119851 ) and in the 
origi nal Baade-Wesselin k method it is neglected 



(see ICaccin et al.lll981 ). However, the Cepheid 
radii estimated by including the AB term (com- 
plete CORS method) in the Eq.((3]) are more ac- 
curate than those based on the origina l Baade- 



2011, and ref- 



Wesselink method (Mo linaro et al 
erences therein). 



4.2. The surface brightness calibration 

To use the complete CORS technique, it is nec- 
essary to calibrate the surface brightness func- 
tion. Here we describe only the main steps of 
the procedure f ollowe d and refer the reader to 
Molinaro et al.l ([201 ll ) and references therein for 
more mathematical details. 

Assuming th e validity of t he q uasi-static 

* Il98fih. for the 



approximation (jOnnembo et al 



Cepheid atmosphere, it is possible to express 
any photometric quantity as a function of ef- 
fective temperature T e and gravity log g. As a 
consequence, considering the surface brightness 
Sy and two generic colors Cy and Cm they 
can be expressed as Sy = Sv(T e ,logg) and 



dj = C^Te, logs), C k i = C kl (T e ,\ogg). If the 
last two equations can be inverted, it is possible 
to express effective temperature and gravity as 
function of the two colors and, consequently, the 
surface brightness becomes: 

S v = SviT^Ci^CulloggiCi^Cu)). (6) 

In general, the invertibility condition is not 
valid over the entire parameter space, since 
the same pair of colors trace different pairs of 
gravity and temperature. However, after an 
appropriate choice of the colors Cy and Cki 
and of their range of variability, it is possible 
to find a local invertibility condition. Using 
grids of models provided by iBessell et al.l (|1998l ) 
( [http://wwwuser.oat.ts.astro.it/castelli/colors/bcp.htmll ) 
and interpolated at the correct metallicity values, 
we succeeded in inverting the previous equations 
for (U-B) and (V-K) colors and for the range 
of parameters typical of Cepheids with period 
P ~ 3 days: i.e. 0.5 < \ogg < 4.5 dex and 
5250 <T e < 6750 K. 

Figure [3] shows the selected region of the theo- 
retical grids in the (V-K) (U-B) color-color plane 
containing the loops of all the analyzed Cepheids. 
For clarity reasons we show only data for a se- 
lected sample of objects. The loops outlined by 
HV12198 and HV12202 are representative of those 
drawn by the majority of Cepheids with the excep- 
tion of We2 and HV12197. The loops of these two 
objects show extreme positions, with the bluest 
and the reddest (U-B) color respectively. This oc- 
currence could be explained invoking the binarity 
of the target or, more likely, the presence of one 
or more blending objects. In this case, depending 
upon the angular separation between the Cepheid 
and the blended companions, relative to the radius 
of the seeing disk, it is possible to oversubtract or 
undersubtract the light of the companion. Simi- 
larly, the loop of HV12203 seems too extended (es- 
pecially along the (U-B) direction) with respect to 
the average of "normal" Cepheids. 

The effect of possible blending on the radius 
calculation has been analyzed in Sec. 15.21 

Within the plotted ranges of colors we have 
been able to fit the relations logT e = logT e (U — 
K,U- B), \ogg = \ogg{V - K, U- B) by means 
of polynomials. The fitted surfaces are shown in 
Figs. |4] and [5] and all the mathematical details 
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are given in Appendix [B] The rms around the fit- 
ted relations amount to 0.00013 and 0.03 dex for 
logT e and logg respectively. 

Using the fitted relations for temperature and 
gravity, we have derived the surface brightness 
from the following equation: 



S v = 42.207- 10.0 log T e - BC 



(7) 



where the constant only depends on the bolo- 
metric absolute magnitude of the Sun, the solar 
constant and on the St efan-Boltzmann constant 
(|Fouaue fc Gierenlll997l) . while BC is the bolomel ,- 
ric correction calculated as a function of the tem- 
perature and gravity through a polynomial fit (see 
Appendix |B]) . Figure |6] shows the fitted surface 
BC = B C (log T e , l ogg) t ogether with the models 
from lBessell et all (|l998f ). 

As a test of our calibration, we have analyzed 
the relation F v — F V (V - K), where Fy = 
4.2207 — 0.1SV is the surface brightness param- 
eter. Figure [7] shows the selected Cepheids plot- 
ted in the Fy-( V-K) plane together w ith the rela- 
tion obtained bv lKervella et all (|2004 using intcr- 
ferometric measurement of 9 Galactic Cepheids. 
The errors on the color (V-K) have been esti- 
mated by considering the scatter around the inter- 
polated light curves, as explained in Sec. 13.11 while 
those on surface brightness have been estimated by 
means of simulations, as described in Appendix [Cl 
below. The plot shows a small discrepancy be- 
tween the fitted relation and the surface brightness 
of NGC 1866 Cepheids, which seems to be system- 
atically bri ghter than the expected values from the 
relation bv iKervella et all (|200i ). Although the 
observed systematic shift is included within the 
uncertainty on red dening, it is importan t to stress 
that the relation bv lKervella et al (2004) has been 
obtained from Galactic Cepheids, which are more 
metal rich than those in NGC 1866. Therefore, 
the observed shift can be due to metallicity differ- 
ences in the sense that less metallic Cepheids are 
more luminous than more metal rich ones. Then 
a metallicity dependence of the surface brightness 
relation is not ruled out. 




log g=4.5\. 



Fig. 3. — Color-color loops of five Cepheids of 
our sample p l otted on the theoretical grids from 
Bessell et all ( 19981) . For clarity only a selected 
sample of objects is plotted (see text for details). 
Locus of constant temperature (dashed lines) and 
of constant gravity (continuous lines) are also plot- 
ted. 



5. Derivation of the radius 

This section contains the procedure followed to 
estimate the radii of the selected Cepheids and a 
comparison with other results present in the liter- 
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Fig. 4. — The surface representing logT e as 
a function of (V-K) and (U-B) obtained from 
the fit of the atmo sphere models (crosses) from 
Bessell et all (|l998l ). 



Fig. 6. — The surface representing the bolometric 
correction BC as a function of logT e and \ogg 
obtained from the fit of the atm osphere models 
(crosses) from lBessell et all (|l998l ). 



log g 




Fig. 5. — The same as Fig. |4]but for logg. 




Fig. 7.— The selected Cepheids of NGC 1866 
(black points) plotted in the plane JFy— (V-K) to- 
gether with the surface br i ghtne ss color relation 
obtained by iKervella et al. | (<2004h (solid line). 
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ature. Here we justify also the value of the pro- 
jection factor p chosen in our analysis. 

5.1. The projection factor 

One of the most uncertain parameter in the 
Baade-Wesselink method is the projection factor 
p, which converts the radial velocity into pulsa- 
tional velocity. It depends both on physical struc- 
ture of stellar atmosphere and the way the radial 
velocity is measured. Furthermore, there is an 
open discussion about its dependence on the pe- 
riod and/or pulsational phase. The question has 
been faced by many authors and the reader can 
find an exhaustive review in Barnes et al. (2009) 
and lStorm et all (|2011al) . 

As mentioned before, the radial velocities 
used in this work have been derived with the 
cross-correlation method. In a recent work, 



iNardetto et al. (l2009h achieved a period depen- 
dent value of the projection factor to correct radial 
velocity obtained by means of cross-correlation. 
Using their relation, namely p = 1.31 — 0.08 logp, 
and the mean value of the period of our sample 
(once first-overtone are fundamentalized accord- 



ing to iFeast fc Catchpold 119971 ). we calculated 
the projection facto r value p = 1.27. T his is the 
same value found by lGroenewege n (2007) by using 
Cepheids with interferometric angular diameters 
and HST parallaxe s. Furthermore i t allow s us to 
be consistent with iMolinaro et al.l (|201lf ). where 



they compare the two values p = 1.36 and p = 1.27 
and found that the latter is the favoured one. 

5.2. Radius calculation 

Using the photometric data and the radial ve- 
locities discussed in Sec. EJ we are able to derive 
the mean radius for each star of our sample by 
using a FORTRAN 77 code. This performs a fit 
of the V magnitude curve, the (V-K), (U-B) color 
curves and the radial velocity curve using a Fourier 
fit, with a number of harmonics fixed interactively 
by the user. Then it solves the CORS Eq. 1 
for the radius at an arbitrary phase both with and 
without the AB term. The mean radius is, finally, 
calculated by integrating the radial velocity curve 
twice. Our results are listed in Tab. |U It con- 
tains the linear radii (in solar units) obtained from 
the CORS method and the angular diameters (in 
mas) calculated from the calibrated surface bright- 



ness Eq.JT]), for the selected sample of Cepheids of 
NGC 1866 analysed in this work. The uncertain- 
ties on both the parameters are also listed and 
are obtained from Monte Carlo simulations as de- 
scribed in Appendix [Cj Note that for HV 12203 
we provide only the radius obtained without the 
AB term (29.5i? Q ). Indeed, the inclusion of this 
term produces an anomalously small radius value 
(25.3i? Q ), as a consequence of the peculiar loop 
of HV12203 in the color-color plane (see Sec. B~2"]) 
that, in turn, generates an unusually large value 
of the AB term (~ 10~ 2 ). The anomalous width 
of the HV 12203 loop has been considered in the 
estimate of the error on the angular diameter for 
this object. 

Table [4] contains also the linear radii obtained 
by other authors for the stars in common with our 
sample. Their values have been rescaled to our 
projection factor value. O ur results are syste mat- 
ically larger than those bylStorm et al. ( 2005T) and 
iGieren. Richtler fc Hilkerl (|1994[ ). whereas they 
are consistent within the uncertainties with those 
by ICote et aL ( 1991 ) , with the exception of HV 
12199 which has a highly undetermined radius in 
his work. 

As noted in Sec. 14.21 the color-color loops of 
We2 and HV12197 are systematically shifted in 
(U-B) with respect to the locus occupied by all 
other Cepheids on the grids of models in the (U- 
B), (V-K) plane. This could be due to the effect 
of overestimating or underestimating the flux in 
blue-ultraviolet bands, due to the difficult pho- 
tometric analysis of crowded regions. To esti- 
mate the possible effect of blending on the de- 
rived radii, i.e. how the presence of a blue unre- 
solved close companion affects the (U-B) color of 
the quoted Cepheids, we decided to shift along the 
(U-B) direction the two extreme loops described 
by We2 and HV12197, to match the region of the 
grids occupied by most of the Cepheids. This 
color shift is of the order of ±0.08-0.09 mag. We 
then re-calculated the radii of these two Cepheids 
with the CORS method, obtaining as a result a 
difference in radii smaller than 1% and 4% for 
We2 and HV12197, respectively. Similarly, for 
the angular diameters the change was 0.3% and 
0.5%, respectively. Therefore, our result is ro- 
bust agaist the effect of contamination of the flux 
of Cepheids by undersubtractcd or ovcrsubtracted 
companions. This is due to the fact that blend- 
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Table 4: Mean linear radii (second column) and angular diameters (third column) of all selected Cepheids 
are listed with their uncertainties estimated from Monte Carlo simulations, as described in the text. The 
linear ra dii, expressed in so lar unit s, listed in the remaining colum ns have been obtai n ed by the following 



linear ra mi, expressed m so lar unit s, listed m tne remaining columns nave been obtaine d by tne 10 
authors: IStorm et al.1 (|2005l) (S05), iGieren. Richtler fc Hilkerl (Il994l) (G94),IC6te et ail (|l991l) (Cm) 



Star 


(R±5R){R Q ) 


(0 ± 56) (/i-arcsec)) 


S05 


G94 


C91 


HV 12197 


31.6 ± 1.7 


5.58 ± 0.08 


23.9±0.6 






HV 12198 


33.1 ± 1.3 


6.00 ± 0.06 


27.5±0.4 


28.3±2.6 


38.3±6.3 


HV 12199 


27.9 ± 1.1 


4.96 ± 0.08 


23.0±1.0 


24.6±3.1 


56.1±19.1 


HV 12202 


29.6 ± 0.5 


5.78 ± 0.08 


26.3±0.9 




25.2±4.4 


HV 12203 


29. 5± 0.7 a 


5.22 ± 0.11 


26.1±1.1 


24.2±4.4 


25.9±4.9 


V4 


30.2 ± 1.3 


5.74 ± 0.04 






29.3±8.1 


V6 


31.8 ± 1.3 


5.10 ± 0.16 








V7 


31.7 ± 2.0 


5.92 ± 0.08 








V8 


31.2 ± 2.3 


5.18 ± 0.06 








We2 


30.7 ± 1.1 


5.68 ± 0.12 








We8 


28.8 ± 1.0 


5.40 ± 0.10 









"Value obtained without the AB term (see text). 



ing affects only the value of the AB term, which 
typically has a small value (see Sec. l4.1|) and influ- 
ence s the radius at most b y about 5% on averag^ll 
(see iMolinaro et al. 201 ll and references therein) . 
In any case, we will take into account this pos- 
sible source of systematic uncertainty, by adding 
it to the random errors on the radius and surface 
brightness (obtained as in Appendix ICj). 

Using the results of our procedure, we investi- 
gated how the selected Cepheids locate with re- 
spect to the Galac tic period-radi u s rela tion re- 
cently obtained bv IMolinaro et ahl ( 2011 ). A vi- 
sual comparison is shown in Fig.[8]and it is evident 
that the Cepheids in NGC 1866 lie on the same 
period-radius relation as Galactic Cepheids. As a 
test we have fitted the equation log R = a + b log P 
to the sample including the 26 Galactic Cepheids 
bv lMolinaro et al" ( 2011 ) and the stars analyzed in 
the present work, with exception of the first over- 
tone pulsators (V6 and V8). The fitting relations 
are the following: 

logi?= (1.120±0.019) + (0.723±0.019)logP (8) 

logi?= (1.110±0.015) + (0.746±0.015)logP (9) 



respectively excluding and including the AB term. 
The inclusion of the Cepheids belonging to NGC 
1866 extends the Period-Radius relation in the 
direction of short periods, giving stronger con- 
straints on the coefficients of the fit. It is impor- 
tant to stress that the equation in the case with 
the AB term remains almost unchanged wit h re- 
spect to that found bv lMolinaro et ah ( 2011 ). but 
with smaller uncertainties on the coefficients of the 
period-radius relation. Concerning the case with- 
out AB term, the fitting e quation is slightly differ - 
ent with respect to that bv lMolinaro et al.1 (|201ll ). 
but in any case consistent within the uncertainties. 

6. Distance to NGC 1866 

Using the values of the angular diameter, de- 
rived from the surface brightness, and the CORS 
Baade-Wesselink linear radius we have estimated 
the distance to NGC 1866 by using the simple 
equation d(kpc) = 2R(UA)/9(mas), where R and 
6 are the mean linear radius and angular diameter 
respectively. 

The typical procedure used in other works 



1 Note that this is not the case for the Ccphcid HV12203 
quoted above whose loop causes a AB term approximately 
1 order of magnitude larger than the typical ones, and in 
turn, a difference of ~15% in the radii. 



(|Fouque fc Gierenll997trStorm et al.l2005Ll2011ah 

consisting in matching the curves of linear radius 
and angular diameter, gives the same result. In 
particular, we have performed a correction for 
eventual phase shift between the two curves and 
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-1 — I — I — I — I — I — I — I — I — I — I — 

AB*0 rms=0.03 

log(R/R„) = (1.110±0.015) + (0.746±0.015)log(P) 




I 1 1 1 1 1 1 1 1— 



AB=0 rms-0.04 

log(R/R o ) = (1.120±0.019) + (0.723±0.019)log(P) 




Fig. 8. — Period-radius relation obtained by com- 
bining the sample o f 26 Galactic Cepheids by 
Molinaro et al. ( 201 lh (empty circles) with those 
selected in the present work (empty triangles). 
The two panels refers to the fit including the 
term AB (upper panel) and excluding it (bot- 
tom panel). The fitting equations and the scatter 
around the fit are also indicated. The first over- 
tone pulsators (filled triangles) are excluded from 
the fit. 



then the distance is calculated as the parameter 
that minimizes the quantity XX -^('AO — ^6*(0i)) 2 , 
where the index i runs over the phases from 0.0 to 
0.8. The phase cut has been performed to avoid 
the influence from possible shocks in the stellar 
atmosphere close to the minimum radius (see e.g. 
Storm et al1l2004h . 



As an example of radial curve matching, we re- 
port the case of We 2. In particular, Fig. |H] shows 
the linear radius of We 2 as a function of the phase 
and the curve of angular diameter corrected for 
phase shift and multiplied for the best distance 
value obtained from the minimization of the func- 
tion cited above. It is evident from the plot that 
the two curves deviate at phases larger than 0.8, 
which have been excluded in the matching proce- 
dure. 

The distances in kpc of the Cepheids analyzed 
in the present work are listed in Tab. [5] The 
corresponding distance moduli, fi = m — M = 
—5 + 5 log[e?(pc)], are also reported in the same ta- 
ble and are plotted in Fig. [TU] as a function of the 
period. 

In our analysis we exclude the two overtone pul- 
sators, whose distances are clearly discrepant (ap- 
parently more distant at a 3cr level) with respect 
to the other considered stars (see Fig. [TU|) . Since 
the Cepheids in NGC 1866 are expected to be all 
nearly at the same distance, this occurrence could 
be due to an imperfect calibration of the surface 
brightness function of overtones throughout the 
grids of model atmosphere. Further analysis is re- 
quired about this point. 

Using a weighted statistic of the 9 remaining 
Cepheids, we obtain a distance of 50.3±0.6 kpc 
corresponding to 18.51±0.03 mag in distance mod- 
ulus and a rms of 0.09 mag. However, if we do 
not consider the weight, the mean distance be- 
comes 51.1±0.6 kpc, equal to a distance modulus 
/i =18.54±0.03 mag with rms of 0.07 mag. Since 
we are confident that our estimate of errors for 
each single Cepheid are correct, we consider the 
weighted result as our best value for the distance 
to NGC 1866. 

The distance obtained in this work is 0.21 
mag l arger than that obtained by IStorm et al. 
(|2005l ). but it is consistent with their recent 
result 18.45±0.04 mag, obtained by calibrat- 
ing the Period-Luminosity relation from the in- 
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frared surface brightness technique, although they 
used a different va lue of the projection factor 
(IStorm et al.ll2011bl) . The distance to NGC 1866 
has been also d erived by usi n g oth er methods. 
As an example, IWalker et al.l (|200ll) used Hub- 
ble Space Telescope V and I photometry of stars 
in NGC 1866 to apply the Main Sequence fit- 
ting technique, obtaining a distance equal to 
18.35±0.05 mag, which is shorter th an our result 



On th e contrary, in a recent work iSalaris et al 



(|2003l ) obtained the distance to NGC 1866 from 
the Red Clump technique. Their analysis gave 
a distance modulus of 18.53±0.07 mag, in agree- 
ment with our result. 

7. Conclusions 

In this work we derived the distance to the 
LMC young populous cluster NGC 1866. To 



technique ( Caccin et al.l 


1981: Rioepi et al.l 1997, 


20001: Ruoddo et al.l 2004 


; Molinaro et al. 201 lh to 



a sample of 11 Cepheids belonging to the clus- 
ter. This method allows to obtain the linear radius 
through the combination of photometric and spec- 
troscopic data and the angular diameter from the 
calibration of the surface brightness function by 
means of grids of model atmosphere. Finally, the 
distance is easily obtained by combining both lin- 
ear and angular diameter in the equation d(kpc) = 
2R(UA)/6(mas). O ur work extends the sample of 



Cepheids studied by IStorm et al.l (|2005l ) thanks to 



new photometric and spectroscopic data. 

We obtained the reddening c orrection by apply- 
ing the technique introduced bv lDean et all (|l978h 



in the (V-I)-(B-V) color-color plane and obtained 
the value E(B-V)=0.06±0.02 mag, which coincides 
with the standard value used for NGC 1866 (see 
Storm et al.| [2005. and references therein). 



As for the projection factor, we used the ap- 
propriate value p=1.27, obtained f rom a recent p- 



factor -log P relation introduced by Nardctto et al 
( 2009 ) and already adopted in [Molinaro et al 



20111) 



The linear radii obtained from the CORS tech- 
nique result to be larger than thos e derived by 



Gicr en. Richtler fc Hil ker (1994) and IStorm et al 



(2005) for sonic stars of our sample. On the con- 
trary, we ar e in good agreem ent with the results 
obtained by ICote et al.l (|199lh . 



We studied also how the selected Cepheids 
place with re spect to the Pe r iod-R adius relation 
obtained by iMolinaro et al.l (|201lh for Galactic 



Cepheids. Our analysis shows that the Cepheids 
in NGC 1866 follow the same linear relation than 
the Galactic ones. 

A weighted statistical analysis gives a distance 
modulus of NGC 1866 equal to 18.51±0.03 mag 
and a rms of 0.09 mag. 

Finally, the distance obtained in this work is in 
agreement with the converging value of 18.50 mag 
for LMC obtained by many authors and described 
in the introduction. 

We are grateful to our anonymous Referee for 
his constructive criticism that helped us to signif- 
icantly improve the paper. Financial support for 
this study was provided by PRIN-INAF 2008 (P.I.: 
M. Marconi). 
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Table 5: The distances of the selected Cepheids expressed in kpc (second column) and their distance modulus 
(third column) obtained from the CORS Baade-Wesselink linear radius and the angular diameters derived 
by the surface brightness. 



Star (d ± Sd) (kpc) Distance modulus (mag) 



HV 12197 


53 


2 


± 


3 





18 


63±0 


12 


HV 12198 


■52 


2 


± 


2 





18 


59±0 


08 


HV 12199 


53 


1 


± 


2 


3 


18 


62±0 


10 


HV 12202 


48 


1 


± 


1 





18 


41±0 


05 


HV 12203 


52 


6 


± 


1 


7 


18 


61±0 


07 


V4 


49 


5 


± 


2 


2 


18 


47±0 


10 


V6 


58 


5 


± 


2 


9 


18 


83±0 


11 


V7 


50 


3 


± 


3 


2 


18 


51±0 


14 


V8 


5G 


3 


± 


4 


2 


18 


75±0 


17 


We2 


51 





± 


2 





18 


54±0 


09 


We8 


50 


2 


± 


2 


1 


18 


50±0 


09 




pha: 



Fig. 9. — Example of matching between linear 
radius curve and angular diameter curve for the 
Cepheid We 2. 



P (days) 



Fig. 10. — Distance moduli (p), with error- 
bars, of fundamental Cepheids (empty circles) 
and first overtones (full points) are plotted as 
function of the period (P) expressed in days. 
Our best distance estimate (solid line) is plot- 
ted together with the r esults by other au - 
thors (point - dashed lines): IStorm et al.l (12011b ). 
Storm et all (l2005h fS05l. IWalker et all (l200lh 
fWOlV iGieren. Richtler fc Hilkerl (I1994T) (G94), 
Cote et al.1 (|l99ll) fC9lV ~ 
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A. Spectroscopic data 

The new spectroscopic measurements for the three Cepheids HV 12197, HV 12199 and We2, were secured 
with FLAMES VLT under the program 074.D-0305. The observations were performed in the GIRAFFE 
mode and employing three different setups, namely HR 11 (with a coverage between 5597 and 5840 A), 
HR 12 (5821-6146 A) and HR 13 (6126-6405 A), with spectral resolutions ranging from 19000 and 24000. 
Five exposures were secured for HR 11, four for HR 12 and three for HR 13 (all the exposures are of 3600 s 
each). 

Each exposure was analyzed independently. The spectra were reduced using the girBLDRS pipeline 
developed at the Geneva Observatorj^J including bias subtraction, flat -fielding, wavelength calibration with 
a reference Th-Ar lamp and the final extraction of each spectrum. The accuracy of the zero-point in the 
wavelength calibration was checked by measuring the position of so me sky-emission l i nes, c ompared with 
their rest-frame position taken from the sky-emission lines atlas by lOsterbrock et al. I (Il996h . No relevant 



shift was detected in the emission lines position, confirming the reliability of the adopted wavelength solution. 
The girBLDRS pipeline performs, in its last step, the measurement of the radial velocity through the 



classical cross-correlation technique ( Tonrv fc Davis 19791) . Heliocentric corrections were applied to each 



radial velocity. The typical uncertainty in the derived velocity of each exposure is of about 0.8 km/s. 

The measurement of t he iron content o f the t hree cepheids was performed following the same approach 



already employed in the iMucciarelli et al.l (|201 lh concerning the chemical analysis of a sample of stars in 
NGC1866. In each spectrum we identify a number of reliable, unblended iron lines (typically 10-15 in each 
spectrum), deriving their abundance through the ^-minimization of the observed line profile with a grid of 
synthetic spectra computed with differe nt Fe abundan ces. The computation of the synthetic spectra were 
performed by means the SYNTHE code ( Kurucz 1993) coupled with the ATLAS9 model atmospheres. 



B. The fit of theoretical grids. 

Here we give the explicit mathematical expression of the surfaces obtained by interpolating logT e and 
log g as a function of the two colors (U-B) and (V-K), and that describing the bolometric correction BC as 
a function of logT e and log<?: 

logT e = ax + a 2 {V - K) + a 3 {V - Kf + a 4 (V — K)(U — B) + (Bl) 
+a 5 (V - Kf{U - B) + a 6 (V - K)(U - Bf + 
a 7 (U - Bf + a s (V - Kf + a 9 (U - Bf + a w (U - B) 



log<? = h + b 2 (V - K) 2 + b 3 (V - K)(U -B) + bi(V - K) 2 (U - B) + (B2) 
+b 5 (V - K)(U - B) 2 + b 6 {U - Bf + b 7 (V - Kf + b 8 (U - Bf + b 9 (U - B) 



BC = a+ c 2 logT e + c 3 (logT e ) 2 + c 4 (log T e ) (log g ) + c 5 (logT e f '(log g) + (B3) 
+c 6 (logT e )(log 5 ) 2 + c 7 (log 3 ) 3 + c 8 (log 3 ) 2 + c 9 logg (B4) 

The coefficients Oj, bi, c% of the previous equations are listed in Tab. |6l The rms of the previous relations 
are 0.00013 dex, 0.0007 dex and 0.0012 respectively. 



1 http: / /girbldrs.sourceforge.net/ 
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Table 6: Coefficients of the polynomial equations obtained from the interpolation of temperature, gravity and bolomctric correction. 



oo 





a 2 


a 3 


a 4 


a 5 


ag 


a 7 


a 8 


ag 


a 10 


3.9858±0.0019 


-0.208±0.005 


0.055±0.004 


0.110±0.008 


-0.031±0.004 


0.057±0.008 


-0.029±0.007 


-0.0093±0.0011 


-0.071±0.007 


-0.121±0.004 


<>1 


f>2 


i>3 


f>4 


65 


!>6 


67 


!>8 


(■9 


b 10 


5.51±0-09 


-4.61±0.1S 


25.9±1.6 


-14.0±0.S 


! 24.5±1.5 


-13.2±1.3 


2.65±0-10 


-21.7±1.4 


-19.4±0.S 




c l 


<=2 


c 3 


<=4 


<=5 


c 6 


c 7 


c 8 


<=9 




-221.0±7.4 


114.9±4.0 


-14.9±0.5 


-21.4±1.8 


2.9±0.2 


-0.106±0.009 


0.00160±0. 00016 


0.39±0.03 


40.2±3.4 





C. Error estimate on radii through Monte Carlo simulations 



To investigate how the errors in the various parameters, involved in the CORS method, influence the 
linear radius estimation, we have performed Monte Carlo simulations. In particular we have considered 
the uncertainties on the colors (V-K) and (U-B), on the reddening E(B-V) and the error in the phase 
matching between the photometric and radial velocity curves. We have run 1000 simulations for each of 
the cited parameters. They consist in varying each analyzed parameter using random numbers extracted 
from a Gaussian distribution with an rms equal to the error on the parameter itself. For all the extracted 
displacements the radius is recalculated and the rms of the obtained values is assumed as the uncertainty in 
the radius due to the error on the analyzed parameter. As for the colors (V-K) and (U-B) we have derived 
their errors from the scatter around the interpolated light curves. Here we do not give all the derived errors 
in the colors, but remind only that typical values are ~0.03 mag for both the (V-K) and (U-B), except 
for few cases with errors of ^0.06 mag because of less sampled light curve in the K and/or U bands. For 
the three Cepheids HV12197, HV12203, and We2, an additional systematic error has been considered for 
the (U-B) color, leading to a total uncertainty of the order of 0.085-0.10 mag. The effect of the reddening 
on the radius estimation has been analyzed by considering the error of 0.02 mag obtained as described in 
Sec. |3l while for the phase matching we have chosen a typical error of 0.01 in phase. This value has been 
chosen to be conservative, because the error on the phase matching, due to the uncertainty on the period 
propagated between the epoch of photometry and that of spectroscopy, is typically sm aller. Furthermore , 
we have taken into account the fact that the star HV 12202 is a spectroscopic binary (jWelch et al.lll99ll ) 
and consequently its photometry and radial velocity can be influenced by the flux and the motion of the 
companion respectively. To be conservative we have doubled the errors on the photometry and the radial 
velocity data. The total effect of the various errors on the radius is always less than 5%, with exception of 
V7 and V8 for which it is -7% (Tab. EJ). 

Monte Carlo simulations have been used also to estimate the errors on the surface brightness. As it 
depends on the two colors (V-K) and (U-B) through the equations given in Appendix [B] we have run 1000 
simulations for each color following the same procedure described before. The derived error on the surface 
brightness has been used to calculate the uncertainty on the angular diameter, obtained from the Eq.Q and 
reported in Tab. |U 
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